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Cosmological evolution of the FRII source population 
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ABSTRACT 

By combining a model for the evolution of the radio luminosity of an individual source 
with the radio luminosity function, we perform a multi-dimensional Monte-Carlo sim- 
ulation to investigate the cosmological evolution of the Fanaroff-Riley Class II radio 
galaxy population by generating large artificial samples. The properties of FRII sources 
are required to evolve with redshift in the artificial samples to fit the observations. 
Either the maximum jet age or the maximum density of the jet environment or both 
evolve with redshift. We also study the distribution of FRII source properties as a 
function of redshift. From currently available data we can not constrain the shape 
of the distribution of environment density or age, but jet power is found to follow a 
power-law distribution with an exponent of approximately -2. This power-law slope 
does not change with redshift out to z = 0.6. We also find the distribution of the 
pressure in the lobes of FRII sources to evolve with redshift up to z ~ 1.2. 

Key words: galaxies: jets - galaxies: evolution - galaxies: intergalactic medium - 
galaxies: luminosity function, mass function - method: statistical 



1 INTRODUCTION 

The formation and evolution of massive galaxies are impor- 
tant components in understanding how the universe evolved. 
It is not clear why AGN formation takes place at the centres 
of some galaxies, but it is believed that this occurs as part of 
the formation of all massive galaxies and jet activity is often 
associated with AGN. Some of the most distant and power- 
ful active galaxies have been observed in the radio band and 
radio astronomy has been closely connected with cosmology 
ever since. 

Radio galaxies with extended lobes are believed 
to influe nce a significant f ractio n of the matte r -filled 
universe dRawlings fc Jarvid |2004 iKronberg et all l200ll : 
iGopal-Krishna fc Wiitall200lf ). From observations in X-rays, 
we know there are significant interactions between radio 
sources and the surrounding gas on scales of several tens to 
hundreds of kpc from the c entral AGN fe.g. lBohringer et al.l 
1 19951 : [Reynolds et alj|2005h - A denser gas environme nt may 
imply a more massive galaxy ijO'Sullivan et al1l200lh . while 
more massive galaxies also contain more massive black holes 
at their centres which can ge nerate more powerful jets (e.g. 
iKormendv fc Richstonef T995h Therefore we can expect that 
there is a connection between the jet properties and their en- 
vironments. The aim of our work is to study the behaviour 
of the properties of the radio jets as a function of their en- 
vironments and of cosmological redshift. 
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iFanaroff fc Rilevl (|l974l ) split extragalactic radio 
sources into two classes depending on their large-scale mor- 
phologies: FRI and FRII. FRI objects have bright cores 
and edge-darkened lobes while FRII objects are edge- 
brightened and contain hotspots. FRIIs are more luminous 
than FRIs with a transition luminosity a round Pi78MHz ~ 
10 25 WHz-'sr" 1 . lOwen fc Ledlowl (|l994h found that the 
transition in radio luminosity extends into the optical and 
FRIIs are preferentially associated with more optically lu- 
minous galaxies. However, the break luminosity between the 
FR classes increases with increasing optical luminosity of the 
host galaxies. The physics leading to the FRI-FRII transi- 
tion is not well understood, but may depend on the gas den- 
sity i n the environment (e.g. lAlexanderlbOOOl ; iKaiser fc Best! 
120071 . hereafter KB07). In our work, we only concentrate 
on the FRII-type objects, because we have good analytical 
models for the evolution of individual objects. 

Various published samples of radio sources 
observed at low radio fre quenc ies, such as 
3CRR llLaing Riley fc Longairl Il983j h 6CE 
dRawlings. Eales fc Lacvf ~ 1200 if ) and 7CRS 

|McGilchrist et al.1 |l99Ch . that are complete down to 
their flux-limits. These samples indicate that the comoving 
density of radio galaxies was higher during the quasar 
era around r edshift z = 2 a s com p ared to the present 
epoch (e.g. IJackson fc Wafll Il999l : IWillott et ait l200ll : 



2004 ). The star and galaxy formation rate 



iGrimes et al.1 

was also considerably higher in the quasar era. The radio 
luminosity function (RLF) has been developed as a tool 
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to describe the densities of radio sources. Willott et al.l 
l|200ll . hereafter W01) generated an RLF at 151 MHz which 
shows a peak at around z — 2. Optical and hard X-ray 
obse rvations of power ful AGN also reveal a similar trend 
(e.g. lUeda et al.|[2003r i. 

Many analytical models for extragalactic classical dou- 
ble radio sources have been published which character- 
ize radio sources in terms of their dynamics and radio 
luminosity evolution . In t his paper we use the model of 
iKaiser fc Alexander] (Il997l . hereafter KA97), which showed 
that the lobes inflated by the jets emerging from the AGN 
grow in a self-simila r fashion. The radio lobe luminosity 
evolut ion is given bv lKaiser. Dennett- Thorpe fc Alexanderl 
(1 19971 . hereafter KDA). The radio synchrotron emission 
of the lobes is due to relativistic electrons spiralling in 
the magnetic field of the lobe. The model of KDA self- 
consistently takes into account the energy losses of these 
electrons due to the adiabatic expansion of the lobes, 
synchrotron radiation and inverse Compton scattering of 
cosmic microwave background pho tons off the electrons. 



Blun dell. Rawlings fc Willottl (|l999l . hereafter BRW) and 
Manolakou fc Kirkl l|2002l . hereafter MK) essentially follow 
the KDA prescription, but differ in the way the relativis- 
tic particles are injected from the jet into the lobe and in 
the treatment of loss terms and particle transport. Thus the 
radio luminosity evolution from these three models show sig- 
nificant differences. In our work, we will generate artificial 
samples by using these three models and compare the results 
with each other. 

One of the most important tools to study the evolu- 
tion of ra d io sou rces is the P-D diagram introduced by 
IShklovskiil (|l963T l, which uses the two fundamental prop- 
erties o f radio sources: their radio luminosity, P, and linear 
size. D. iBaldwinl (|l982l) pointed out that it is analogous to 
the Hertzsprung-Russell diagram for stars. The source dis- 
tribution in the P-D plane is a result of the intrinsic evo- 
lution of individual sources and the cosmological evolution 
of the source population as a whole. Since we may assume 
that source lifetimes are considerably shorter than the age of 
the universe, the P-D diagram has been used to pl ace some 
constraints on the evolu tion of individual sources (I Baldwin! 
1 19821 ; IN eeser et alii 19951 ) and to look for consistency between 
data and models (KDA; BRW; MK). In our work, we will 
investigate the evolution of radio sources across cosmolog- 
ical epochs, therefore we take into account redshift at the 
same time and extend to a 3-dimensional P-D-z data cube. 

The purpose of this work is to use the evolutionary 
model of individual radio sources together with the RLF 
to generate artificial samples with large numbers of sources. 
From these artificial samples we can find the best fitting 
parameters describing the radio sources and their environ- 
ments, how the jet properties are distributed and how they 
evolve over c osmological time scal e s. Ou r approach differs 
from that of IKaiser fc Alexander] l| 19991 . hereafter KA99) 
who assume a 'birth function' to describe the probability 
of the radio source progenitors becoming active and turn- 
ing into radio sources. Their intrinsic luminosity evolution 
is determined by the properties of their jets and the environ- 
ments the progenitors are located in at some cosmological 
epoch. We extend this work by directly using the RLF from 
W0 1 as a description of the birt h function. 

iBlundell fc Rawlingsl i|l999l) also investigate the trends 



of radio galaxy properties with redshift. The main difference 
between their work and ours is that they use a model for 
the radio hotspots leading to steeper tracks in the P-D dia- 
gram. They also assume a 'birth func tion' to investig a te the 
distribution of the w hole population . Barai fc Wiital (120061 . 
hereafter BW06) and lBarai fc Wiital (|2007l . hereafter BW07) 
test the same three evolutionary models for FRII sources 
we use here. They show that none of them fit the observa- 
tional data, but again they only take into account the 'birth 
function' instead of the RLF. The use of a birth function 
in KA99, BRW, BW06 and BW07 is designed to provide a 
good fit of the respective model to the high-luminosity end 
of the RLF. Our approach uses the RLF directly without the 
need for a 'good guess' at the birth f unction. We compare 
our re sults with those of BW06 and IBlundell fc Rawlingsl 
(|1999T ). 

In section 2 we describe the RLF generated by W01. In 
section 3 we describe the observational samples we use in 
our work. We then describe the P-D-z data cube we use to 
compare the samples in section 4. Key features of the KA97 
and KDA models are discussed in section 5 and we sum- 
marize our Monte-Carlo method in section 6. In section 7 
we introduce the 3-dimensional KS test which we use in the 
sample comparison. We show our simulation results in sec- 
tion 8 and compare the results from different radio emission 
models in section 9. A discussion of the simulation results 
is performed in section 10 and finally we give our main con- 
clusions in section 11. Throughout the paper we are using a 
cosmological model with Ho = 71kms _1 Mpc -1 , S7m = 0.3, 
Q A = 0.7. 



2 RADIO LUMINOSITY FUNCTION AT 
151 MHZ 

To construct our artificial samples we need to know the rel- 
ative number of objects with a given radio luminosity at 
a given redshift. This information is provided by the radio 
luminosity function (RLF) denned as the number of radio 
sources per unit co-moving volume and per unit logarithm to 
base ten of luminosity at a given redshift, p(P, z). Several de- 
terminations of p(P, z) at various observing frequencies are 
available in the literature. The analytical model we use to 
connect source parameters with observable source properties 
models the radio emission of the lobes without the hotspots 
at the end of the jets. The emission from the hotspots is most 
important at high frequencies and may even dominate the 
total emission. To minimize the effect of the hotspot emis- 
sion we use the RLF at 151 MHz c ompiled by W01 on the 
basis of 356 so urces from the 3CRR (|Laing. Riley fc Longairl 
Il983ft. 6CE llRawlings. Eales fc Lacvl l200ll ) and 7CRS 
l|McGilchrist et al.lll990r i flux-limited samples. 

W01 model the RLF as the sum of two distinct pop- 
ulations allowing independent redshift evolution. The low- 
luminosity population, pi, contains a mixture of FRI-type 
sources and the lowest luminosity FRII-type objects. The 
high luminosity population, ph, contains only FRII-type 
sources. The total RLF is then p(P, z) = pi + ph- 

The low-luminosity population is modelled as a 
Schechter function, 



Pi = fi{z)pw ip^ 



' leXP (£. 



(1) 
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Figure 1. The adopted radio luminosity function corresponding 
to model C of W01 for Q M = 0.3, Q A = 0.7 and Q k = 0. Dashed 
and dotted lines show the low-luminosity and high-luminosity 
population respectively. The solid lines show the sum of both 
components. 



where pio is a normalization term. At luminosities P below 
the break luminosity Pi* the RLF approximates a power-law 
with slope —a\. The low-luminosity population decreases ex- 
ponentially above the break. The normalization of pi evolves 
with redshift through 



Mz) = (l + z) k 



(2) 



up to a maximum redshift Zio beyond which f\ remains 
constant. Here we use model C of W01 and so we adopt 
log pio = -7.523, on = 0.586, log Pi* = 26.48, fci = 3.48 and 
zio = 0.710. 

The high-luminosity population is parameterized in a 
similar way as 



Ph = fh(z)pho ( -p— 



exp 



(3) 



where the exponential cut-off is now located below the 
break luminosity Ph* and the power-law with slope — ah ex- 
tends above the break. The number of sources in the high- 
luminosity population is modelled as rising up to z = Zho 
and then decreasing at higher redshifts as 



h = exp 



/h 



exp 




for 2 < Zho 



for 2 ^ 2 h0 . 



(4) 



The relevant constants are logpho — —6.757, = 2.42, 
log P h * = 27.39, 2 h0 = 2.03, z hl = 0.568 and z h2 = 0.956. 

In this paper we only model sources of type FRII. W01 
do not distinguish between the FR classes in their determi- 
nation of the RLF. While ph contains only FRII-type ob- 
jects, the exact composition of pi in terms of FR class is 
not known. Here, we simply assume that 40% of the sources 
contributing to the low luminosity part of the RLF are of 
type FRII. We find below that this assumption allows for a 
good fit of the properties of our artificial samples to those 
of the observed samples. However, the fraction of FRII-type 



sources in pi may be a function of redshift and/or lumi- 
nosity. In fact, in Section 8.4.2 below we show that the 
observed sample with the lowest flux limit we use in this 
paper, the 7CRS sample, requires the FRII fraction to in- 
crease with increasing redshift. The advantage of using the 
RLF directly rather than a birth function is that we can 
include the weaker FRIIs in the low-luminosity population 
while the birth functions used by BRW, BW06 and BW07 
only consider the high-luminosity population which exclu- 
sively contains FRIIs. 

W01 compute the RLF for a cosmological model with 
H = SOkms^Mpc -1 , M = 0, = and fi k = 1. 
We adopt the cosmological parameters consistent with the 
WMAP results, H = 71 kms" 1 Mpc" 1 , Q M = 0.3, fi A = 
0.7 and fik = 0. Hence we need to convert the R LF, p, to 
the c orrect cosmological model using the relation (|Peacockl 
fl985h : 



Pi (A, 2)— = p 2 (P2,2) 



(5) 



2 d2 

where P is the luminosity derived for a measured flux and 
redshift 2 in a specific cosmological model, while V is the 
comoving volume. The indices refer to two different cosmo- 
logical models. 

The comovi ng distance in a given cosmological model is 
fe.g.lHogdli99^ 

dz' 



D M {z) = 



Ho J 



y/Sl M (l + Z'Y + fi k (l + 2') 2 + «A 



(6) 



The measured flux of a source at redshift 2 corresponds to 
different luminosities in different cosmological models and 
they are related by 



The comoving volume is 
dV = 4ttD 2 m AD m . 



(8) 



We use the above relations to translate the RLF of W01 into 
our adopted cosmological model. We show the adopted RLF 
in Figure [T] 



3 COMPLETE SAMPLES OF RADIO 
SOURCES 

Flux-limited samples contain all radio sources within a well- 
defined sky area with radio fluxes above a specified limit at 
the observing frequency. The term complete sample refers 
to a flux-limited sample of radio sources in which the cos- 
mological redshifts of the host galaxies and the linear size 
of the lobes of all sample members are also measured. Fi- 
nally, since the model used here only describes sources with 
FRII-type morphologies, we also need to know the classi- 
fication of each sample member in terms of FR class. W01 
used three complete samples, 3CRR, 6CE and 7CRS, for the 
determination of the RLF. We use the same samples here. 
However, in addition to the measured redshifts and radio 
luminosities of all sample members necessary for the con- 
struction of the RLF, we will also use their measured linear 
sizes. The himinosity-redshift distribution of the FRII-type 
objects contained in the three samples is shown in Figure 
[2] Note the range in radio luminosity covered by the sam- 
ples at a given redshift. Below we will also use the complete 
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Figure 2. The radio luminosity-redshift plane for the 3CRR, 
6CE and 7CRS samples. The different symbols identify sources 
from different samples: 3CRR (pluses), 6CE (asterisks) and 7CRS 
(squares). 



BRL sample l|Best. Rottgering fe Lehnertll 19991 ) to estimate 
the goodness of fit of our models. We do not directly use the 
BRL sample to constrain our models because its flux limit is 
defined at 408 MHz, significantly higher than the observing 
frequency of the RLF we use to construct our models. Below 
we briefly describe the samples we use in this paper. 

The Third Cambridge Re vised Revised sample (3CRR; 
lLaing. Riley fe Longairl [1983) covers the largest sky area 
(4.23 sr) of all the samples used here, but it is also the 
shallowest. Its flux-limit of Si 78 > 10.9 Jy is defined at 
178 MHz rather than 151 MHz, the observing frequency of 
the RLF we use. The sample contains 146 FRII-type sources. 
To convert the flux limit as well as the fluxes of indi- 
vidual objects from 178 MHz to 151 MHz W01 simply as- 
sume a typical spectral index of -0.8 and we follow this 
prescription. The flux limit therefore translates to 12.4 Jy 
at 151 MHz. We use the redshifts, angular sizes and lu- 
minosities of sample members listed online by Willott at 

www.science.uottawa.ca/~cwillott/3crr/3crr.html. 

The Sixth Cambridge sample (6CE) defined by lEalesI 

i|l985l ) is b ased on the 6C survey. I t was re-selected and 
updated bv lRawlings. Eales fc Lacvl (|200lh . The radio flux 
at 151 MHz of sources in this sample falls in the range 
2.0 sC Sibi < 3.93 Jy. Note that 6CE is the only sample 
that has an upper as well as a lower flux limit. The area 
covered by the sample is 0.103 sr. There are 52 FRII-type 
sources in this sample. We use the r edshifts and luminosities 
of sam ple members as published in lRawlings. Eales fc Lacvl 
J200lh while the angular sizes of the sources are taken from 
lEalesf l|l985h . 

The Seventh Cambridge Redshift Survey (7CRS) is a 
combinatio n of the sub-divisions I ,H and III of the original 
7C survey (|McGilchrist et alJll99Ch . The sample is observed 
at 151 MHz and contains every source with a flux-density 
S151 ^ 0.5 Jy in three regions of the sky covering a total of 
0.022 sr. The data relevant for this w ork for the third fi eld of 
the sample can be found in Table 9 o f lLacv et all <|l999l ). The 
data for the first and se cond fields of the sa mple is referred 
to, but not published, bv lGrimes et al.l (|2004T ) and can be ob- 



tained from www-astro.physics.ox.ac.uk/~sr/grimes.html. 
The classification of sources in the 7C-I and 7C-II sub-fields 
of the sample into FR types has not been published. 

We obtained the relevant data from the VLA archive 
to repeat the FR classification. The raw data was reduced 
using standard routines within AIPS. Within the 7C-I and 
7C-II fields we find 11 sources of type FRI, 52 sources with 
an FRII-type morphology and 11 compact sources. The com- 
pact sources are unresolved in all of the archival data down 
to an angular size of 0.5 arcsec. To prevent problems in the 
construction of our artificial samples, we exclude any artifi- 
cial sources with linear sizes corresponding to this smallest 
resolvable angular size or smaller. We list our source classi- 
fica tion for the 7C-I and 7C-II fi e lds in Table 1 and 2. 

iBest. Rottgering fc Lehnertl (|l99 9) define a complete 
sample (BRL) at an observing frequency of 408 MHz and 
a flux limit of 5 Jy. For a typical spectral index of -0.8 this 
flux limit at 408 MHz translates to a limit close to that of the 
3CRR sample at 178 MHz. Thus the BRL sample allows us 
to estimate the similarity of samples drawn from the same 
parent population with similar selection criteria, but in dif- 
ferent parts of the sky. However, given the large difference 
in observing frequency we do not attempt to fit the BLR 
sample directly with our model. The BRL sample contains 
178 sources, 124 of which have an FRII-type morphology. 



4 THE P-D-Z DATA CUBE 

We use the RLF derived from observations to constrain the 
number of sources in our artificial samples with a given radio 
luminosity at a specified redshift. Therefore, by construc- 
tion, our artificial samples agree with the two-dimensional 
distribution of radio luminosities, P, and redshifts, z, the 
RLF, or P-z diagram, of the observed complete samples. 
However, the complete samples also contain the measured 
linear sizes of the radio lobes and therefore represent a three- 
dimensional data cube with the axes radio luminosity, P, lin- 
ear size, D, and redshift, z. To agree with the observational 
constraints, our artificial samples must provide a good fit 
to the observed source distribution in the full P-D-z data 
cube. 

Previous work has concentrated on two-dimensional 
projections of the full data cube of varying subsets of the 
three samples we discuss here (e.g. Subramanian & Swarup 
1990; Neeser et al. 1995; KA99; BRW; Arshakian & Longair 
2000; W01; MK; Grimes et al. 2004; BW06 and BW07). In 
this paper we use the full three-dimensional data cube to 
constrain the models to avoid the loss of information asso- 
ciated with the projection process. 

BRW, BW06 and BW07 also considered the radio spec- 
tral index as a fourth parameter. They found that the ob- 
served distribution of the spectral index is not well fitted. 
They use the same models for the evolution of radio sources 
that we employ in this paper. The models significantly re- 
strict the possible range of the spectral index and this effect 
naturally explains the difficulties encountered by BW06 and 
BW07. Here we do not attempt to reproduce the observed 
distribution of spectral indices as this would either result in 
the same problem or require the introduction of an addi- 
tional model parameter. We will return to this point below. 
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Table 1. Data for the 7C-I sub-field. Column 1 is the source name, column 2 is the redshift, column 3 
is the linear size in kpc, column 4 is the logarithm of the radio luminosity at 151 MHz and column 5 is 
the morphology classification with 1 indicating FRI, 2 indicating FRII. c indicates an unresolved, compact 
object. Column 6 is the program code of the VLA observations used for the FR classification. 



Name 


z 


Linear size (kpc) 


logioCPisi/WHz-^sr- 1 ) 


Morphology 


Program code 


5C6.5 


1.038 


186.4 


26.31 


2 


AB371 


5C6.8 


1.213 


47.4 


26.83 


c 


AR335 


5C6.17 


1.05 


396.1 


26.56 


2 


AB371 


5C6.19 


0.799 


68.3 


26.51 


2 


AR335 


5C6.24 


1.073 


11.1 


26.68 


2 


AB766 


5C6.25 


0.706 


198.2 


26.08 


1 


AL355 


5C6.29 


0.72 


92.4 


25.94 


2 


AL355 


5C6.33 


1.496 


124 


26.65 


2 


AB371 


5C6.34 


2.118 


66.4 


27.13 


2 


AL355 


5C6.39 


1.437 


214.5 


26.59 


1 


AL355 


5C6.43 


0.775 


31.6 


26.16 


1 


AB371 


5C6.62 


1.45 


271 


26.93 


1 


AL355 


5C6.63 


0.465 


370.6 


25.66 


2 


AR477 


5C6.75 


0.775 


112 


25.94 


2 


AR365 


5C6.83 


1.8 


119 


27.2 


2 


AR365 


5C6.78 


0.263 


1459.7 


25.62 


2 


AB667 


5C6.95 


2.877 


113 


27.55 


2 


AR335 


5C6.160 


1.624 


53.5 


26.88 


2 


AB371 


5C6.201 


0.595 


76.9 


26.12 


1 


AR365 


5C6.214 


0.595 


216.4 


25.98 


2 


AR477 


5C6.217 


1.41 


103.1 


27.13 


2 


AR335 


5C6.233 


0.56 


48.3 


26.07 


2 


AL355 


5C6.237 


1.62 


23.8 


27.11 


c 


AR365 


5C6.239 


0.805 


616.4 


26.2 


2 


AB766 


5C6.242 


1.9 


42.6 


27.06 


2 


AR365 


5C6.251 


1.665 


50.7 


26.7 


2 


AR365 


5C6.258 


0.752 


2.4 


26 


c 


AR365 


5C6.264 


0.831 


40.6 


26.28 


2 


AL355 


5C6.267 


0.357 


23.7 


25.16 


2 


AB766 


5C6.279 


0.473 


183.6 


25.55 


1 


AR365 


5C6.282 


2.195 


8 


27.04 


c 


AR365 


7C0221+3417 


0.852 


140.8 


26.8 


2 


AL355 


5C6.286 


1.339 


140.3 


26.65 


2 


AL355 


5C6.288 


2.982 


7.3 


27.6 


c 


AR335 


5C6.287 


2.296 


103.9 


27.57 


2 


AR335 


5C6.291 


2.91 


4.4 


27.57 


c 


AB667 


5C6.292 


1.241 


41.4 


26.72 


2 


AB371 



5 MODELS FOR THE DYNAMIC AND 

EMISSION EVOLUTION OF INDIVIDUAL 
SOURCES 

The large-scale structure of FRII-type sources is inflated by 
powerful jets accelerated in the vicinity of the supermassive 
black holes at the centre of the AGN inside the host galaxies. 
The jets end in strong shocks which accelerate electrons to 
relativistic velocities and may increase the strength of the 
magnetic field. The magnetized plasma passing through the 
jet shock subsequently inflates the lobe or cocoon . The basic 
dynamical picture was fi rst proposed bv lScheuei] (|l974l ) and 
iBlandford fc Reesl l|l974h . 

iFallel (|l99lf ) and KA97 showed that if the jet is in 
pressure-equilibrium with its own lobe, then the expansion 
of the lobe and the bow shock in front of it is self-similar. 
Here we use the model of KA97 to describe the dynamics 
of the lobes. We summarize the most important features of 
this model below. 

The growth of the jet length is essentially determined 
by a balance of the ram pressures of the jet material and 



that of the medium surrounding the host galaxy which is 
pushed aside by the jet. X-ray observations of the hot gas 
in the vicinity of elliptical galaxies, in galaxy groups and 
in galaxy clus ters find density distr ibutions well fitted by 
/3-models (e.g. iFukazawa et al]|2004T ). In KA97 the density 
distribution outside the core radius ao is approximated by a 
power-law: 

where d is the radial distance from the AGN at the centre 
of the density distribution and po is the density at the core 
radius. The exponent j3 is constrained by observations to the 
range < (3 < 2. 

We assume that the gas in the vicinity of the AGN 
has a non-relativistic equation of state, V x = 5/3, while the 
lobes only contain magnetic fields and relativistic particles, 
r c = 4/3. For a jet providing a constant power Qq for a time 
t and inflating a lobe with an aspect ratio Rt we then find 
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Table 2. Same as Table 1, but for the 7C-II sub-field. 



Name 


z 


Linear size (kpc) 


logioCPiSi/WHz^sr- 1 ) 


Morphology 


Program code 


5C7.7 


0.435 


13.4 


25.58 


1 


AL355 


5C7.8 


0.673 


320.3 


26.36 


2 


AB371 


5C7.9 


0.233 


440.8 


25.37 


2 


AB371 


5C7.10 


2.185 


170.3 


27.54 


2 


AB371 


5C7.15 


2.433 


16.4 


27.36 


c 


AR335 


5C7.17 


0.936 


691.5 


26.21 


2 


AL0401 


5C7.23 


1.098 


235.2 


26.63 


2 


AB371 


5C7.25 


0.671 


6.3 


25.78 


c 


AB667 


5C7.47 


1.7 


1.7 


26.79 


c 


AB371 


5C7.57 


1.622 


634.7 


26.79 


2 


AB371 


5C7.70 


2.617 


13.7 


27.75 


2 


AR365 


5C7.78 


1.151 


187.6 


26.99 


2 


AR365 


5C7.79 


0.608 


1863.9 


25.77 


2 


AL355 


5C7.82 


0.918 


358.3 


26.28 


2 


AL355 


5C7.85 


0.995 


227.4 


26.64 


2 


AL0401 


5C7.87 


1.764 


94.8 


27.17 


2 


AR335 


5C7.95 


1.203 


486.8 


26.65 


2 


AL0401 


5C7.106 


0.264 


104.8 


25.28 


1 


AB371 


5C7.111 


0.628 


80.4 


26.29 


1 


AB371 


5C7.118 


0.527 


76.3 


26.06 


2 


AL355 


5C7.125 


0.801 


120.2 


26.15 


2 


AB371 


5C7.145 


0.343 


93.2 


25.31 


2 


AB371 


5C7.170 


0.268 


97.1 


25.19 


2 


AB371 


5C7.178 


0.246 


121.6 


25.15 


1 


AB371 


5C7.194 


1.738 


16.8 


27.3 


2 


AB371 


5C7.195 


2.034 


22.1 


27.12 


2 


AB371 


5C7.205 


0.71 


107.3 


26.34 


2 


AB371 


5C7.208 


2 


146.7 


27.27 


2 


AL0401 


5C7.223 


2.087 


42.2 


27.07 


2 


AL355 


5C7.242 


0.992 


389.9 


26.22 


c 


AL355 


5C7.245 


1.61 


100.2 


27.23 


2 


AB371 


5C7.269 


2.218 


61.6 


27.21 


2 


AL355 


5C7.271 


2.224 


9.6 


27.06 


2 


AR365 


5C7.400 


1.883 


491.9 


27.14 


2 


AR365 


5C7.403 


2.315 


15.8 


26.96 


c 


AR365 


7C0825+2446 


0.243 


375.8 


24.94 


1 


AR365 


7C0825+2443 


0.086 


122.1 


24.86 


2 


AB371 



which we assume to initially follow a power-law with expo- 
nent p between a low and high energy cut-off represented 
by the Lorentz factors of the least and most energetic elec- 
trons, 7 m in and 7max, respectively. To simplify the equations 
and the model we choose p = 2 for all sources. We also set 
7min = 1 and 7max = 10 5 . Note that the exact value of the 
power-law limits does not significantly influence the model 
results as long as 7 m i n -C 7max- A summary of all relevant 
model equations can be found in KB07. 

Our choice of p limits the range of spectral indices pre- 
dicted by the model to the narrow range from —0.5 to — 1. 
For young sources, synchrotron energy losses are most im- 
portant, while for old objects inverse Compton losses dom- 
inate and both processes lead to a steep spectrum. For 
'middle-age' sources between these extremes, the relativistic 
electrons are not much affected by radiative energy losses 
and so the spectrum flattens (KB07). In objects with ob- 
served spectral indices not conforming to this pattern a dif- 
ferent acceleration regime may be applicable and lead to a 
different value of p. Our choice of a fixed value for p im- 
plies that we cannot expect the spectral index distribution 
of our artificial source populations to agree with that of the 
observed complete samples. However, the introduction of an 
unknown distribution for p into the model would consider- 



the lobe length, D, and the pressure inside the lobe, p c 
\Poao / 



27c? p o, o 3 .-4-fl 



16(5 - /3) : 



R^( Po a^t^Qt r , (10) 



where 



_ . 64i%(5 -f3f , 



^81tt(8-/3 + 3J?| 

Note here that the model of the source dynamics only de- 
pends on the combination poa^, but not on pa and ao sepa- 
rately. For convenience we therefore introduce the parameter 
A = poOo- in our model we initially assume (3 — 2 as this 
considerably simplifies the model equations. We will discuss 
the model with /3 = 1.5 later. 

KDA97 extended the dynamics model of KA97 to in- 
clude the calculation of the radio emission expected from the 
lobe. The model takes into account the energy losses of the 
relativistic electrons due to the adiabatic expansion of the 
lobe, synchrotron radiation and inverse Compton losses from 
scattering of the cosmic background radiation. These losses 
modify the energy distribution of the relativistic electrons 
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Figure 3. The flow diagram of the Monte-Carlo process we use 
in this paper. 



ably complicate the analysis of our results and introduce 
another model parameter. 

The combined model of KA97 and KDA in the form 
specified here depends on five model parameters, z, Rt, A, t 
and Qo. A complete set of these parameters fully determines 
the linear size of the lobes, D, and their radio luminosity, P. 
Therefore we can take the observed distribution of sources 
in the P-D-z space to constrain the underlying distribution 
of the density of their environments described as A, their 
age t and their jet power Qo- In the following we explain 
the practical implementation of this process as well as our 
assumption for the distribution of lobe aspect ratios, Rt- 



6 MONTE-CARLO SIMULATION 

The RLF, p(P,z), gives the comoving number density of 
radio sources at a cosmological redshift z with a given ra- 
dio luminosity P. The number of sources within the ranges 
z to z + dz and logP to logP + dlogP is given by n = 
p(P,z) dlogPdV, where dV is the comoving volume sam- 
pled between z and z + dz. The relevant formulae are pre- 
sented in Section 2. We consider redshifts from z — 4 to 
z = with a step size of dz — 0.004 and radio luminosities 
at an observing frequency of 151 MHz from log P — 24.5 to 
log P = 29 in steps of d log P = 0.01. These ranges cover the 
entire source population in all three observed samples we 
consider. We only calculate n for combinations of P and z 
for which the corresponding observable flux exceeds the flux 
limit of the observed complete sample we use to compare 
our artificial sample with. 

After calculating n at all positions defined by our ranges 
and step sizes, we normalize n such that its maximum is 



equal to unity. We can now interpret n (P, z) as the dis- 
tribution of the probability to find a source in an artificial 
sample with a given combination of redshift and radio lu- 
minosity. To construct our artificial sample we then iterate 
through all allowed combinations of P and z and generate 
a random number, r, uniformly distributed between and 
1. If r ^ n(P,z) for the chosen P and z, then a source 
with this radio luminosity and redshift is included in the 
artificial sample. We can increase the sample size by repeat- 
ing the same process a fixed number of times at all allowed 
combinations of P and z. 

Once a source with a given combination of radio lumi- 
nosity and redshift is included in the artificial sample, we 
need to calculate the linear size of its lobe to determine its 
position in the P-D-z data cube for comparison with the 
observations. For this, we need the aspect ratio of the lobe, 
Rt, the age of the source, t, and the parameter describ- 
ing the external density distribution, A. Fixing the values 
of these three parameters uniquely determines the length of 
the lobes, D, and the jet power, Qo. In other words, fixing 
the values of three parameters uniquely determines the val- 
ues of the remaining two. For each source in the artificial 
sample we will choose random values for Rt, t and A from 
distributions of these parameters discussed in the following 
in order to assign values to their linear sizes and jet power. 

It is difficult to accurately determine the aspect ratio 
of the lobe of a radio source. The lobe should ideally be de- 
tected all the way back to the central AGN which suggests 
the use of a low observing frequency to minimize the effects 
of radiative energy losses. At the same time we need suf- 
ficient angular resolution to resolve the lobe perpendicular 
to the jet axis. This is best achieved at high observing fre- 
quencies. In practice a balance must be found between these 
requirements and so an observationally determined distribu- 
tion of Rt is not available in the literature. For simplicity we 
choose a uniform distribution across the range 1.3 ^ Rt ^ 6, 
where the limits are motivated by observed values (Leahy & 
Williams, 1984). Below we discuss the minor effects on our 
results of changing this distribution. 

If we assume that all radio sources have a common max- 
imum lifetime t max , then the distribution of their ages is 
given by a uniform distribution extending from t = to t = 
tmax- Spectral ages of some objects reach a few 10 s yea r s (e.g . 
I Alexander fc Leahvl [l987h . but recently iBird et all (|2008f) 
found an average age of around 10 7 years for a sample of 
FRII-type objects located in low redshift galaxy groups, sug- 
gesting a maximum lifetime of a few 10 7 years. The maxi- 
mum lifetime depends on the availability of fuel for the jet- 
producing AGN and it would therefore be surprising, if t max 
was the same for all so urces. However, th e spread around 
the average age found bv lBird et alj l|2008h is comparatively 
small and for simplicity we assume that i max is indeed the 
same for all sources at a given redshift. We will show below 
that t m ax defined in this way may be a function of redshift. 

The parameter A describing the density distribution 
in the environment of a source can, in principle, be de- 
termined from X-ray observations. FRII-type objects seem 
to be preferentially located in isolated galaxies or galaxy 
groups. With P — 2, we expect A b etween 10 17 kgm~ 1 an d 
10 18 kgm _1 for individu al galaxies (Fukaz awa et alJl2004l ). 
For group environments iJetha at al. ( 20071 ) find values for 



A in the range 2 x 10 18 kgm _1 and 4 x 10 19 kgm _i . Note 
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that all these determinations predict somewhat high val- 
ues for A, because we are using f3 = 2 at the moment. For 
/3 = 1.5 the values of A are more cons istent with the va l- 
ues given bv lFukazawa et al.l l|2004 ) and lJetha at al.1 (|2007h . 
The gas haloes of many galaxies and galaxy groups show 
flatter slopes implying smaller values of A. Again there is no 
determination of the complete distribution of A for the envi- 
ronments of radio sources available in the literature. Hence 
we draw random values for A from a uniform distribution 
extending from A = to A = A max , where the maximum 
value is to be determined from the models and may be a 
function of redshift (see below). 

For every allowed combination of the radio luminosity 
and the redshift plus the randomly chosen values for _Rt, t 
and A we can now proceed to calculate a linear size D of 
the lobe and jet power Qo for each source to be included in 
the artificial sample. The parameter distributions used here 
allow for a wide range of possible values for Qo- Not all val- 
ues for Qo are acceptable. Sources with weak jets, i.e. small 
jet powers, will develop turbulent, ra ther than laminar jets 
and th is will result in FRI-type lobes. iRawlings fc Saunders! 
(1991) suggest that the transition between the FR types oc- 
curs close to Qo = 10 37 W based on radio observations of 
lobes and optical line emission from the AGN itself. A simi- 
lar lower limit for the jet power of sources with an FRII-typc 
morphology was derived by KB07 using the model we em- 
ploy here as well. An upper limit on Qo may be given by 
the Eddington luminosity of the most massive black holes. 
In our model we require that 10 37 W < Qo < 5 X 10 40 W. 
Model parameters that imply a jet power outside this range 
are rejected. For rejected sources we generate new sets of 
model parameters, _Rt, t and A, in the way described above 
until an acceptable value for Qo is found. 

Note here that the mathematical form of the distribu- 
tions for the model parameters only indirectly influences the 
distributions of these parameters in the final artificial sam- 
ple. Combinations of model parameters must lead to ac- 
ceptable results in the sense that not only must the result- 
ing source have the correct radio luminosity and redshift, 
its jet power must also fall within the specified range. The 
resulting parameter distributions will in general not follow 
the distributions they are originally drawn from, because 
many possible parameter combinations will be rejected as 
described above. In large enough artificial samples the choice 
of uniform distributions over other mathematical functions 
therefore only weakly influences the final distributions of the 
model parameters. 

Once an acceptable combination of model parameters 
is found, we can calculate the corresponding length of the 
lobes, D. This is the physical length of the lobe and in order 
to compare to the observed complete samples we need to 
take into account that the lobes may be projected as their 
main axis is oriented at an angle 9 to our line of sight. The 
observed lobe length is therefore given as D Q b = Dsin^, 
where the random orientation 8 is distributed as sinOdO 
between 8 = and 8 = n/2. Before including a source in 
the P-D-z data cube of our artificial sample, we choose a 
random orientation angle and 'project' the lobes into the sky 
plane. For simplicity, hereafter we denote the D OD as just D. 

The process of generating our artificial samples de- 
scribed above is summarized in the flow diagram in Figure 
[3] Our artificial samples contain of order 6000 sources. This 



is a much larger number than is contained in the observed 
samples. However, we found that such a large number is re- 
quired to arrive at reasonably smooth source distributions 
in the P-D-z data cubes and also to minimize the influence 
of the initial model parameter distributions on the final dis- 
tributions. 



7 KOLMOGOROV-SMIRNOV TEST 

Our artificial samples define a source distribution within the 
P-D-z data cube equivalent to the observed samples. We use 
the three-dimensional version of the Kolmogorov-Smirnov 
(KS) test to compare our model results with the observa- 
tions. 

The classical one-dimensional Kolmogorov-Smirnov test 
makes use of the probability distribution of the quantity 
Dks, defined as the largest absolute difference between the 
one-dimensional cumulative distributions of two samples, 
where one or b oth samples can be continuous or discrete. 
iPeacockl <|l983h extends this idea to a two-dimensional test 
by making use of the maximum absolute difference between 
two distributions, when all four possible ways to cumulate 
data following the directions of the coordinate axes are con- 
sidered. For the comparison of a sample with n members 
with a continuous distribution Peacock's test requires that 
the cumulative distributions should be calculated in all 4n 2 
quadrants of the plane defined by, 

{x <X i>y < Yj), (x <X t ,y> Yj), (x > X h y < Yj), 

(x>X i ,y>Yj), (12) 

for all possible combinations of the indexes i and j from 1 
to n. Here, Xi and Yi denote the coordinates of individual 
me mbers of the discrete sample. 

iFasano fc Franceschinil (|l987T ) show that it is sufficient 
to consider only the four quadrants defined by each indi- 
vidual data point in the discrete sample, i.e. i = j. This 
reduces the total number of quadrants in which the distri- 
butions are accumulated to 4n. Furthermore, the extension 
of this methodology to three dimensions is straightforward, 
provided all eight quadrants defined by each individual data 
point are considered in deriving the largest difference, -Dks, 
of the cumulative distributions. 

Using the three-dimensional KS test we can now com- 
pare the fit of various artificial source distributions within 
the P-D-z data cube arising from different models with that 
of the observed samples. However, the method does not as- 
sign a formal goodness-of-fit measure because the statistics 
of -Dks in the three-dimensional KS test is in general not 
known. Also, the selection criteria for the observed sam- 
ples, in particular their flux limit, prevent the population 
of certain parts of the P-D-z cube. To assess which model 
provides an acceptable fit to the observations, we separately 
construct the statistics of -Dks for each model calculation. 
For this we generate a large number of sources for a specific 
model as detailed above. We then repeatedly draw random 
subsamples from these model sources with a total source 
number equal to that of the observed comparison sample. 
We calculate the total difference Dks of the cumulative dis- 
tributions of the subsample and the parent model sample. 
In this way we build up the probability distribution of -Dks 
for this particular model. Based on this statistic we assign a 
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probability P(-Dks) to the value Dks calculated for the ob- 
served sample. P(Dks) is then the probability that the ob- 
served sample arises as a subsample drawn randomly from 
all radio sources in the universe, if that population is de- 
scribed by the model. While this technique may guide us in 
the selection of models, it clearly cannot identify the sta- 
tistically most likely model because the distribution of Dks 
will be different for each individual model. 



8 RESULTS 

We now have all the ingredients to generate our artificial 
samples and compare them to the observed samples. In do- 
ing so our general approach is to start with the simplest 
models and only modify these as necessary to achieve a bet- 
ter agreement with the observations. It is not feasible to 
show the full three-dimensional source distributions used in 
the comparisons. In order to present our results, we therefore 
plot two-dimensional projections of the P-D-z cube. Sources 
in observed samples are presented as individual crosses while 
the large numbers of sources in the artificial samples are 
shown as density contours in these plots. The contours en- 
close areas of 1%, 10%, 40% and 80% of the maximum den- 
sity in each plot. 

Note that in the following we always show the radio 
luminosity density of sources in the 3CRR sample and as- 
sociated artificial samples at 178 MHz rather than 151 MHz. 
While we use the RLF at 151 MHz in our models to calcu- 
late relative source numbers as detailed above, we calculate 
the luminosity for our artificial sources at 178 MHz and use 
the appropriate flux limit for this frequency when comparing 
with the 3CRR sample. 

8.1 Model A 

The simplest model we can build within the framework de- 
scribed above has uniform distributions with fixed upper 
limits for the source age and the density parameter. We as- 
sume a fixed £ max = 2.5 x 10 7 years and investigate a range 
of possible A max . None of the A max tried leads to a probabil- 
ity significantly different from zero as measured by the KS 
test. Varying t max instead of A max leads to the same result. 

As an example, in Figure [4] we show a comparison of 
the resulting artificial samples with A max = 1 X 10 18 kgm -1 
with the observed samples. Clearly high luminosity sources 
located at high redshifts are too large in our artificial sam- 
ples compared to the observed sources. The sources in the 
artificial samples show a trend of increasing size with in- 
creasing luminosity until the trend is reversed at the highest 
luminosities. The radio luminosity of a source in the KDA 
model is mainly governed by the jet power, Qo. However, 
a larger value of Qo also implies faster growth of the lobes 
and, all other parameters being equal on average, it is not 
surprising that we find a trend of radio luminosity with size 
in our artificial sample. 

The reversal of this trend at the hi ghest luminosities 
has been noted b y many authors (e.g. lOort et al.l Il987l ; 
iNeeser et af] [l995) for observed samples. In the artificial 
samples it is caused by the flux limits of the samples in com- 
bination with the decreasing luminosity of older, and there- 
fore larger, sources. Clearly the evolutionary model used 



here predicts that this effect alone is not sufficient to ex- 
plain the observed trend reversal and some other effects must 
contribute to reconcile the model with the observations (see 
also K A99) . We will investigate such effects in the next two 
sections. 

Modified models with a steeper luminosity evolution of 
individual sources do not need to invoke such additional ef- 
fects (BRW and MK). However, they may be less consistent 
with the overall source distribution (BW06) . We shall apply 
these models in the same way to the data as the combined 
models of KA97 and KDA in section 9. 



8.2 Model B 

In this and the following sections we investigate which ad- 
ditional factors may lead to a better fit of the predicted 
distribution of sources in the P-D-z cube with the obser- 
vational data. We have shown above that in our modeling 
framework we need an additional effect to explain the appar- 
ent shortening of the lobes of sources with the highest radio 
luminosities located at the highest redshifts. In this section 
we introduce a redshift dependence for the upper limit of 
the density parameter A. In the following section we do the 
same for the upper limit of the source age distribution. 

We fix the maximum source age to i max = 5 x 10 7 years 
and introduce a variable upper limit for the distribution of 
A such that A max (2:) = A max (0) (1 +zy. The best agree- 
ment between artificial and observed samples is achieved 
for 3CRR with A max (0) = 3.7 ± 0.2 x lO^kgrn" 1 and 
if) = 5.9 ± 0.2, giving a probability that the observed sam- 
ple is drawn from a population described by the model 
of P(D KS ) = 48%. The values for 6CE are A max (0) = 
3.8t°oi x 10 18 kgm-\ i/> = 5.8 ± 0.2 and P{D KS ) = 56%. 
The limits on the model parameters give their value where 
P(Dks) halves compared to its maximum. The fit of the 
model to the 7CRS sample is much worse and the values 
for P(-Dks) are very low. This is related to problems with 
the used RLF for this sample. We discuss this issue in Sec- 
tion 8.4.2. In any case, the maximum of P(-Dks) for 7CRS 
occurs close to A max (0) = 3.7 x 10 18 kgm -1 and ip = 6, con- 
sistent with the results for the other two samples. A model 
with A max (0) = 3.7 x 10 18 kgm _1 and ip = 5.8 provides 
an acceptable fit to both 3CRR (P(Dks) = 40%) and 6CE 
(P(Dks) = 41%). The comparison of this model and the 
observed samples is shown in Figure [S] 

The much better fit of Model B compared to Model 
A is explained by the considerably higher average density 
of the source environments at high redshifts. The density 
parameter A has only a moderate influence on the radio 
luminosity of a model source compared to the jet power. 
However, a high value of A efficiently reduces the expansion 
speed of the lobes. Hence the strong evolution of A max with 
z in Model B ensures that sources at high redshift remain 
smaller for longer than their low redshift counterparts. This 
allows a better fit of the high luminosity / high redshift part 
of the source population. 

8.3 Model C 

Instead of invoking denser average source environments at 
higher redshifts to reduce the average size of the highest 
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Figure 4. Projections of the P-D-z data cube for observed samples with source density contours from Model A. The parameters we use 
here are i max = 2.5 X 10 7 yr and A max = 1 X 10 18 kgm _1 . Crosses indicate FRII sources in the observed samples. Gray scales indicate 
the number density of the artificial samples. 



luminosity sources, we can also reduce the average source 
lifetime for increasing redshift. For this we fix A max (,z) at 
6.4 x 10 17 kgm -1 for all z and introduce a variable maximum 
age as t max (z) = t max (0) (1 + z) . The best agreement is 
found for t max (0) = 2.T^% x 10 7 yr and cj> = 2.4 ± 0.2 for 
3CRR with P(Dks) = 76%. The 6CE sample gives an al- 
most identical result with i max (0) = 2.7^t 'g x 10 7 yr and 
4> = 2.4 ± 0.1 where P(Dks) = 59%. Again, the compari- 
son with 7CRS does not produce good fits for the reasons 
we will discuss in Section 8.4.2, but the maximum of the 
probability occurs close to the parameter values for the two 
other samples. The maximum source age at low redshifts is 
reassuringly close to the value of t he average source lifetime 
of 1.2 x 10 7 yr recently derived bv lBird et all (|2008t ) for an 
independent sample. Figure [6] compares the model results 
with the observed data. 

We assume a uniform distribution of A in Model C. 
The real distribution of A is unknown and it does not 
have to be uniform. Instead of an uniform distribution, 
we also introduce a Gaussian distribution with its peak at 
3.2 x 10 17 kgm" 1 and a width of a = 1 x 10 17 kgm -1 . The 
distribution is truncated so that there are no values of A 
below or above 6.4 x 10 17 kgm _1 . This distribution cov- 
ers the same range of A, but is of course different from an 



uniform distribution. The agreement with the same £ max (0) 
and <j> are P(D KS ) = 77% for 3CRR and P(D KS ) = 52% 
for 6CE. For the extreme assumption of a fixed value of 
A = 3.2 x 10 17 kgm _1 for all sources, we can still find a 
comparable agreement with P(Dks) = 73% for 3CRR and 
P(Dks) ~ 55% for 6CE. Clearly, we can not rule out these 
distributions compared to an uniform distribution. We will 
discuss this point in detail in section 10. 

Both models B and C can provide adequate fits to the 
observed samples 3CRR and 6CE. Given the results detailed 
above, we cannot formally decide whether it is a reduced 
maximum lifetime or a denser environment that limits the 
lobe sizes of sources with the greatest radio luminosities at 
the highest redshifts accessible to the observed samples. Of 
course, a combination of the two effects is also not ruled out. 
However, the result shows that with currently available sam- 
ples it is not necessary to introduce yet more complicated 
models to explain the distribution of sources in the P-D-z 
data cube. We will concentrate on Model C in the discus- 
sion of parameter dependencies, alternative models and the 
further implications of our results. This does not indicate 
that we prefer Model C over Model B, but is simply done to 
avoid confusion. 
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Figure 5. Projections of the P-D-z data cube comparing observed samples with the source density contours resulting from Model B. 
The parameters we use here are t max = 5 X 10 7 yr, A max (0) = 3.7 X 10 18 kgm — 1 and tp = 5.8. Crosses indicate FRII sources in the 
observed samples. Gray scales indicate the number density of the artificial samples. 



8.4 Potential problems 

8.4-1 The power-law exponent, (3 

(3 is the power-law exponent which indicates how fast the 
environment density decreases away from the center AGN. 
In previous sections, we use (3 = 2, but most studies of the 
environments of low redshift radio galaxies in X-ray band 
imply that (3 could be closer to 1.5 than 2. Here we apply 
(3 = 1.5 in the KDA model to check how much it will affect 
our results. 

With a smaller value of (3, the density of the envi- 
ronment decreases more slowly and the speed of growth 
of the lobe is slower. In this case, it is not surprising to 
have a higher maximum jet age. The best fitting parame- 
ters we find are i max (0) = 1.1 x 10 s yr, (f> = 2.6 for fixed 
A max = 3.5 x 10 8 kgm -1 ' 5 . The godness of fit decreases 
slightly to P(Dks) = 72% for 3CRR and P(£>ks) = 41% for 
6CE. The agreement is still good while the value of t max (0) is 
still consistent with the range discussed in Section 6. There- 
fore, the value of (3 does not significantly change our simu- 
lation results and conclusions. We will still concentrate on 
the models with (3 = 2 in the following sections. 



8-4-2 The axial ratio, Rt 

In our construction of the artificial samples we have used a 
fixed uniform distribution for the aspect ratio of the lobes 
Rt- Since Rt affects both the size and the radio luminosity 
of the lobes, it is reasonable to ask how much our results 
depend on our assumptions for Rt- 

If we replace the uniform distribution of Rt over the 
range 1.3 to 6 we used so far with a constant Rt for all model 
sources, then the distribution of lobe sizes within our arti- 
ficial samples becomes slightly narrower. The average size 
also shifts somewhat depending on the value assumed for 
Rt- This effect is illustrated in Figure [7] There is no notice- 
able restriction of the luminosity range. For Model C with 
the best fitting parameters, but a constant value for Rt, 
we find that the probability P(Dks) decreases to 34% for 
3CRR (Rt = 3.4) and remains constant at 59% for 6CE 
(Rt = 2.8). We conclude that changing the range or distri- 
bution of Rt will only have a minor effect on our results. 
This is particularly encouraging because the alternative evo- 
lutionary models for individual sources we introduce below 
assume fixed values for Rt. 
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Figure 6. Projections of the P-D-z cube for the observed samples with source density contours from Model C. The parameters we use 
here are i max (0) = 2.7 X 10 7 yr, A max = 6.4 X lO^kgm" 1 and <j> = 2.4. Crosses indicate FRII sources in the observed samples. Gray 
scales indicate the number density of the artificial samples. 



8.4.3 The RLF 

The adoption of an RLF allows us to include the low- 
luminosity FRII sources in the artificial sample, and is a 
more direct method to constrain the source distribution than 
the introduction of a birth function. However, in general a 
RLF requires more parameters than a simple 'birth func- 
tion'. The changes in these parameters may also change the 
simulation results. Here we u se the general i zed lu minosity 
function(GLF) constructed by I Grimes et all (|2004h for com- 
parision with the RLF of W01 used so far. 

iGrimes et all i|2004h consider both radio luminosity and 
optical luminosity of the AGN and introduce a parameter 
a encoding the L151 — iom correlation and a parameter 
j3' encoding scatter about this correlation. The GLF based 
on a and /3 1 can generate a smoother RLF than the RLF of 
W01. However the total source counts predicted by the GLF 
do not change very much compared to the RLF of W01 as 
IGrimes et al.l (|2004i ) show in their Figure 11. If we substitute 
the GLF into Model C with the same best fitting parameters 
and a modified FRII fraction of 50% in the low-a population, 
we can still find a reasonable agreement with P(-Dks) = 34% 
for 3CRR and P(D K s) = 26% for 6CE. The use of an RLF 
introduces many parameters, and the changes in these pa- 



rameters do affect the fitting result. However within a rea- 
sonable range, we can still get good agreement between ar- 
tificial and observed samples with a modified RLF. Mean- 
while, the conclusion that the FRII properties evolve with 
redshift do not change with the choice of RLF. 

The FRII fraction in the low-luminosity population is 
another parameter introduced in the RLF that can affect 
the fitting result. The value of 40% comes from the simu- 
lation. If we apply a smaller fraction, for example 20%, we 
will find too few low-luminosity sources. The agreement for 
3CRR and 6CE drops below 1%. Therefore, although some 
fainter surveys at low frequency such as the Bologna surveys 
indicate the FRII fraction at the low luminosities could be 
smaller, we still need to choose 40% as the appropriate value 
for our work for fitting the 3CRR, 6CE and 7CRS samples. 

8.4.4 The 7CRS sample 

In all our models discussed above we noted that our ap- 
proach cannot provide an adequate fit to the source distri- 
bution of the 7CRS sample. The main problem appears to 
be the different relative number of sources at low redshifts 
and high redshifts in the various samples. This may be due 
to a changing mix of sources with different FR morphology 
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Figure 8. Projections of the P-D-z cube of the observed and artificial samples with a modified low-luminosity population of the RLF. 
The artificial sample is generated by Model C with the best-fitting parameters. 



below the break in the RLF. So far we have assumed that a 
constant fraction of 40% of the sources contributing to the 
low luminosity part of the RLF of W01 are of type FRII, 
irrespective of redshift. A closer look at the three observed 
samples reveals that this assumption may be too simplistic. 

Within a range from 10 24 WHz -1 sr" 1 to 
10 26 WHz-'sr 1 at 151 MHz the 3CRR sample con- 
tains 20 sources in a redshift range 0.005 ^ z ^ 0.16 of 
which 12 show an FRII-type morphology. In the same lumi- 
nosity range, 7CRS contains a total of 33 sources spanning 
a redshift range from z = 0.086 to z = 0.775 of which 26 are 
of type FRII. These numbers suggest an increased fraction 
of FRII-type sources within the low luminosity population 
at higher redshifts. The 6CE sample does not fit into this 
trend with 3 FRII-type sources out of a total of 7 within 
the luminosity range and at redshifts intermediate between 
3CRR and 7CRS. However, the number of sources in this 
sample is very small. 

Based on these numbers we test whether a redshift de- 
pendence of the FR mix helps to improve the model fit to 
the 7CRS data. We set the FRII fraction within the low lu- 
minosity population of the RLF of W01 equal to 0.3 + 2 up 
to z — 0.7 and keep it at unity at higher redshifts. With 
the modification of the RLF the agreement of Model C with 
the 7CRS data increases, but not significantly. The value of 



P(Dks) increases for 3CRR to 84%, while it decreases for 
6CE to 47%. A comparison of the model with the observa- 
tional data is shown in Figure [8] 

Changing the FR mix in the low luminosity part of the 
RLF does not have the desired effect of improving the fit 
with the 7CRS data. However, the presently available data 
do not allow us to decide whether a change with redshift 
in the composition of the radio source population is taking 
place or not. 

From Figure [8] we see that there are too many sources 
in our artificial sample around z = 2 and Pi 51 = 
10 27 WHz -1 sr _1 compared to the 7CRS sample. In fact, 
this problem, that the RLF of W01 overpredicts the num- 
ber of sources with this radio luminosity at this redshift 
compared to 7CRS, is discussed by W01 (see their Figures 7 
and 8). The 7CRS data is consistent with no further evolu- 
tion of the high luminosity part of the RLF beyond red- 
shift 2=1. Implementing such a constant RLF at red- 
shifts beyond z = 1 in our model dramatically improves 
the agreement between our Model C and the 7CRS data to 
P(Dks) = 17%. However, this RLF reduces the agreement 
between Model C and the 3CRR sample to P(D K s) = 0.1% 
while P(-Dks) for the 6CE sample drops below 0.1%. There- 
fore it is unlikely that this modification of the RLF is nec- 
essary. The 7CRS sample covers a very small sky area com- 
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Figure 9. Projections of the P-D-z cube of the observed and artificial samples with both the low and high luminosity populations of 
the RLF modified. The artificial sample is generated by Model C with the best fitting parameters. 



pared to the 3CRR and 6CE samples. This may imply that 
the members of this sample are not fully representative of 
the entire source population. Comparisons of the artificial 
samples arising from this modification of the RLF with the 
observations are shown in Figure [9] 



8.5 Comparing the 3CRR sample with an 
equivalent sample 

The BRL sample, as discussed in section 3, is roughly equiv- 
alent to 3CRR with a similar flux limit and source number. 
Thus the BRL sample and the 3CRR sample can be con- 
sidered as two samples drawn from the same universe with 
similar selection criteria, but from different areas of the sky. 
The comparison between these two samples indicates the 
agreement between two samples drawn from the same par- 
ent population and can therefore give us some idea of how 
well our models agree with the observations. 

The BRL sample is observed at 408 MHz, so we use a 
constant spectral index of —0.8 to convert the radio luminos- 
ity of BRL sources to 178 MHz. We ignore the sources whose 
luminosity is below the flux limit given by 3CRR. Compar- 
ing the 3CRR sample and the BRL sample in the P-D-z 
data cube, we get Dks = 0.2662 from the 3-D KS test. If 



we choose an artificial BRL sample with the same source 
number as the real BRL sample from the best-fitting artifi- 
cial 3CRR sample and compare with the real 3CRR sample, 
-Dks is mostly between 0.1 to 0.2. As the smaller values of 
Dks indicate a better fit, this shows that our best-fitting 
artificial sample agrees with the observations to a degree 
similar to the agreement between similar observed samples 
drawn from the same source population. Figure \W\ shows 
the 3CRR sample and the BRL sample in the P-D-z data 
cube. 



9 COMPARISON BETWEEN MODELS OF 
RADIO LOBE EVOLUTION 

Before discussing our results in more detail, we now as- 
sess how much they depend on our particular choice for 
the model of the evolution of individual sources. The model 
describing the dynamics and expansion of the radio lobe, 
KA97, essentially relies on the condition of ram pressure bal- 
ance between the jet material and the receding am bient gas 
in fron t of it. This condition was first introduced bv lScheuerl 
( 1974) and has formed the basis for virtually all subsequent 
models of the dynamical evo lution of radio sources w i th an 
FRII-type morphology (e.g. iBegelman fc Cioffil [l98i iFallej 
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Figure 12. Projections of the P-D-z cube of the observed samples with source density contours generated by the BRW model in the 
framework of Model C. The parameters we use here are t max (0) = 5.0 X 10 7 yr, A ma x = 1.6 X 10 18 kgm - 1 and <f> = 2.5. Crosses indicate 
FRII sources in the observed samples. Gray scales indicate the number density of the artificial samples. 



Il99ll : iNathl Il995l : IChvzvlll997h . While the details of the 
derivation of the lobe dynamics differ, the basic principle 
is the same. As we do not consider in detail the evolution of 
individual objects, we therefore do not use another dynam- 
ical source model. 

Based on the dynamical model of KA97, two alterna- 
tive descriptions for the emission properties of the lobes have 
been formulated by BRW and MK. The main difference be- 
tween these models and KDA is the treatment of the radi- 
ating electrons as they propagate from the ends of the jets, 
or hotspots, into the lobes. The KDA model assumes that 
the energy distribution of the electrons is described by a 
power-law and that they only suffer adiabatic losses during 
the transition from the hotspot into the lobe. BRW considers 
the effect of synchrotron energy losses within the hotspot by 
introducing a doubly broken power-law for the energy distri- 
bution. MK describe in detail the diffusion of the electrons 
from the hotspot into the lobe taking into account radiative 
energy losses and possible re-acceleration. Figure [TT] shows 
a comparison of the luminosity evolution of a single radio 
source as described by the three models. The luminosity 
evolution is steeper for the two alternative models with the 
BRW model the steepest. In the following we will describe 
both models in more detail and then substitute them for 



the KDA model in our formalism. BW06 and BW07 use the 
same models in their comparisons. 

When comparing the results arising from the use of 
the three different evolutionary models, we employ Model 
C with A max = 1.6 x 10 18 kgm -1 , consistent with the origi- 
nal formulation of BRW and MK. However, to allow a direct 
comparison with our previous results we set j3 = 2, 7 m i n = 1 
and 7 m ax = 10 s . Also note that the two models in their orig- 
inal form require a fixed value for Rt- As shown in Section 
8.4.1, this should not lead to major changes of the results. 
Finally, we also want to compare the model predictions with 
the 7CRS data. In order to do so we use the modified RLF as 
described in Section 8.4.2. For the 3CRR and 6CE samples 
we continue to use the unmodified RLF. 

9.1 The BRW model 

The BRW model takes the lobe dynamics from the KA97 
model. This sets the pressure in the lobe, p c , according to 
equation (|10[) . As in KDA, the pressure in the hotspot in 
front of the end of the jet is higher than p c by a constant 
factor. With the assumption of energy equipartition in the 
hotspot, the hotspot pressure fixes the strength of the mag- 
netic field and the normalization of the energy distribution 
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Figure 13. Projections of the P-D-z cube of the observed samples with source density contours generated by the MK model within the 
framework of Model C. The parameters we use here are t max (0) = 4.0 X 10 7 yr, A max = 1.6 X 10 18 kgm - 1 and <f> = 2.5. Crosses indicate 
FRII sources in the observed samples. Gray scales indicate the number density of the artificial samples. 



of the electrons at a given point in time. BRW assume fur- 
ther that the hotspot has a fixed physical size and that it 
takes electrons between 1 yr and 10 5 yr to escape from the 
hotspot into the lobe, independent of any other model pa- 
rameter. These timescales give rise to two breaks in the elec- 
tron energy spectrum. Their positions are calculated using 
the standard equipartition formulae. Below the low energy 
break the electron energy spectrum is modeled as a power- 
law with an exponent of p = 2. Above the high energy break, 
p = 3 is assumed while in between the breaks any value 
2 ^ p ^ 3 can be adopted. In our implementation of the 
BRW model we use p = 2.5 in the intermediate region. Af- 
ter the electrons escape into the lobe, they preserve their 
double power-law energy distribution, but the breaks in the 
power-law evolve with time. 

Figure \H\ shows the result of replacing the KDA model 
in our method with the BRW model. The best agreement of 
Model C is achieved with i max (0) = 5.0 x 10 7 yr and </> = 2.5. 
The artificial sample contains too few large sources at high 
redshift compared to 3CRR. The agreement is very low at 
P(Dks) ~ 0.1%. The steep luminosity evolution predicted 
by the BRW model leads to older and therefore large sources 
dropping quickly below the flux limit of the observed sam- 
ples, in particular at high redshifts. In a sense the BRW 



'over-solves' the problem of the lack of large objects at high 
redshift. This was also noted by MK and BW06. 

The results for the 6CE sample and the 7CRS sam- 
ple, with the modified RLF, are better, since these samples 
contain fewer large sources at high redshifts. However, the 
ratio of small to large sources is still too large at high red- 
shifts. The agreement for 6CE is P(Dks) = 19% and for 
7CRS -P(-Dks) ~ 0.1%, in both cases lower than for the 
KDA model. We conclude that the steep luminosity evolu- 
tion predicted by the BRW model provides a worse fit to the 
observational data than the flatter evolution arising from the 
KDA model. 

9.2 The MK Model 

The MK model is similar to the BRW model as it also as- 
sumes that the jets end in hotspots of a constant physical 
size which also contain a magnetic field of constant strength. 
The radiating electrons are accelerated to an initial power- 
law energy distribution within the hotspot. In contrast to 
the BRW model, MK exactly follow the subsequent evolu- 
tion of the energy distribution under the influence of adi- 
abatic and radiative energy losses as the electrons diffuse 
through the hotspot into the lobe. MK adopt a diffusive 



18 Y. Wang and C. R. Kaiser 




12 3 4 

Redshift, z 



Figure 10. 3CRR and BRL objects in the P-D-z diagram at 
178 MHz. Cross symbols refer to 3CRR while squares refer to 
BRL. We convert the BRL sources from 408 MHz to 178 MHz by 
using a common spectral indexes of 0.8 



transport model in which the mean square distance trav- 
eled by individual electrons on their way from the hotspot 
into the lobe is proportional to t a , where t is the time since 
their acceleration. A value of a = 1 corresponds to the diffu- 
sion regime and we adopt this here, since the sub-diffusion 
regime, a < 1, makes the model comparable to the BRW 
model and the supra-diffusion regime, a > 1, approximates 
the KDA model. 

The MK model takes into account adiabatic and ra- 
diative losses of the electrons during the diffusion process. 
However, MK find that the adiabatic losses lead to luminos- 
ity evolution of their sources in disagreement with observa- 
tions. To avoid this problem we adopt their model B, with all 
corresponding parameter settings, which neglects adiabatic 
losses. We also follow MK in setting the ratio of the diffu- 
sive transport time and the cooling time of an electron, their 
parameter t, to 2 x 10 -3 . After the electrons have diffused 
into the lobe, the MK model describes the further evolution 
of the electron energy distribution in the same way as the 
KDA model. An example for the luminosity evolution using 
the MK model in the way described here is shown in Figure 

El 

Replacing the KDA model with the MK model in our 
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Figure 11. Evolutionary track for three sources with different 
models. The upper curves are for Qo = 1 X 10 40 W and z = 1.5, 
the centre curves are for Qo = 1 x 10 39 W and z = 0.5, and the 
lower curves are for Qo = 1 X 10 38 W and z = 0.2. Each of the 
solid, dashed and dotted curves refer to the tracks given by KDA, 
MK and BRW models respectively. The pluses on the curves are 
time markers denoting source lifetimes of 1,10, 20, ...90 Myr. 

method we find the best fitting Model C with the parame- 
ters i max (0) = 4.0 ± 0.2 x 10 7 yr and <j> = 2.6± 0.1 for 3CRR 
at P(Dks) = 53%. For 6CE we find P{D KS ) = 58% for 
tmax(0) = 4.2 ± 0.3 x 10 7 yr and </> = 2.5 ± 0.1. The values 
of the model parameters almost agree between the two sam- 
ples and Figure [T3] shows a comparison of the artificial and 
observed samples for i max (0) = 4.0 x 10 7 yr and 4> — 2.5 re- 
sulting in P(Dks) = 41% for 3CRR, P(Dks) = 50% for 6CE 
and P(Dks) ~ 1% for 7CRS, with the modified RLF. The 
fitting result for 7CRS is not as good as that for the KDA 
model, because the evolutionary tracks of the MK model 
are steeper and predict more sources with small linear size 
at the high redshift / high luminosity end. 

For the 3CRR and 6CE samples, the MK model pro- 
vides a level of agreement between the artificial and observed 
samples similar to that of the KDA model. The somewhat 
higher value for i max (0) we find for the MK model is par- 
tially due to the large value of A max we use for this model 
compared to the KDA model. However, the steeper lumi- 
nosity evolution of the MK model also requires somewhat 
longer lifetimes of the sources to accommodate the larger 
objects. With the currently available data we cannot decide 
which of the two models provides a better description. We 
continue to focus on the KDA model because its mathemat- 
ical formulation is simpler than that of MK. 



10 DISCUSSION 

We produce artificial samples of radio-loud AGN with an 
FRII-type morphology and compare their properties with 
those of samples of observed objects. The artificial samples 
are consistent with the observed ones provided that there is 
some cosmological evolution of the radio source population. 
Our models require either that the density of the source 
environments increases on average for increasing cosmolog- 
ical redshift, or that the lifetime of the jet flows decreases 
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Figure 14. Final A distribution of the artificial samples gener- 
ated by Model C in two redshift ranges. The upper plane is for 
0^2^ 0.2 and the lower plane is for 1.0 $J z $J 1.2. JVtotal i s 
the total source number in each redshift bin. The solid line is for 
3CRR, the dotted line for 6CE and the dashed line for 7CRS. 
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Figure 15. Final distribution of t in the artificial samples gen- 
erated with Model C in two redshift ranges. The upper plane is 
for 0^2^ 0.2 and the lower plane is for 1.0 $J z $C 1.2. JV to tal i s 
the total source number in each redshift bin. The solid line is for 
3CRR, the dotted line for 6CE and the dashed line for 7CRS. 



with increasing redshift. A combination of both effects may 
also be at work. To simplify the following discussion of the 
properties of the artificial samples, we concentrate on model 
C, implying the cosmological evolution of the jet lifetimes. 
However, we do not consider this model superior to Model 
B. 

The steeper luminosity evolution of the radio lobes de- 
scribed by the model of BRW leads to a worse agreement 
of the model predictions with observations. The model pro- 
posed by MK provides comparable results to those of the 
KDA model. BW06 come to a similar conclusion. BW07 
consider modified models taking into account the increase 
of hot-spot sizes with jet lengths. The modified BRW and 
MK models produce better fits which are at least as good as 
the KDA model, while the modified KDA model produces 
a worse fit. Here we do not attempt to construct a new, 
improved model for individual FRII-type objects and so we 
concentrate on the original KDA model which is consider- 
ably simpler and mathematically less complex. 

In constructing our artificial samples we choose random 
values for the density parameter A and the source age t from 
uniform distributions. The jet power Qo is then adjusted to 
give the model source the correct radio luminosity. Not ev- 
ery possible combination of the set of three parameters A, t 
and Qo is allowed because of restrictions on the magnitude 
of Qo and the selection criteria of the observed sample we 
compare with. It is therefore not clear a priori whether the 
distributions of A and t amongst the objects within the final 
artificial sample are also uniform. Any deviation from a uni- 
form distribution in the final sample may reveal a genuine 
property of the source environments or source ages in the 
universe. 

We first investigate the distribution of A in the arti- 
ficial sample arising from Model C. Figure [14] shows the 
binned distribution of A for two different redshift ranges, 
^ z ^ 0.2 and 1 ^ z ^ 1.2. At low redshifts the distri- 
bution is uniform within the fluctuations arising from the 
finite number of sources in the artificial samples, regardless 



of which observed sample we compare with. At high redshifts 
the distribution of A in the artificial sample mimicking the 
3CRR sample is biased towards large values. The distribu- 
tions within the artificial samples compared with 6CE and 
7CRS remain uniform except for a drop in the bin of the 
smallest density. This behaviour is caused by the flux limits 
of the samples. Objects located in denser environments are 
more luminous. At low redshifts the flux limit of all samples 
corresponds to such low luminosities that the entire radio 
source population with an FRII-type morphology is repre- 
sented in the samples. However, at higher redshifts the flux 
limit excludes sources in less dense environments as their 
luminosity is too low. Samples with a lower flux limit like 
6CE and 7CRS obviously suffer less from this problem. 

The result for the final distribution of A arising from an 
initially uniform distribution implies that our current data 
is insufficient to constrain the distribution of the environ- 
mental properties of FRII-type sources. The uniform dis- 
tribution 'survives' unchanged through the source selection 
process. However, the distribution of A in the universe itself 
may be uniform and thus, by chance, we have selected the 
correct distribution. To test this we also generate artificial 
samples with a Gaussian distribution and also a fixed value 
of A. As discussed in Section 8.3, the artificial samples using 
these two alternatives agree with the observed samples to a 
similar degree as the standard Model C with a uniform dis- 
tribution. The final distribution of A in the artificial sample 
drawing random values from a Gaussian distribution is also 
Gaussian. By construction, all the sources in the final sam- 
ple with a fixed value of A are assigned this one value. From 
these results we conclude that the distribution of A in the 
universe cannot be constrained with the currently available 
data. 

The distribution of source ages in our artificial sam- 
ples for two redshift bins is shown in Figure [151 The distri- 
butions at low redshift are uniform with a slow decrease 
for old ages. The decrease is caused by two effects. For 
Model C the source lifetime is decreasing with increasing 
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redshift. Hence within the respective redshift ranges there is 
also by construction a range of source lifetimes. Sources to- 
wards the high redshift end of the range cannot contribute 
to the bins of t corresponding to the oldest sources. The 
second effect is again due to the flux limit of the sam- 
ples. The radio luminosity of sources decreases as they get 
older. Therefore older so urces are less likely to be included 
in the samples (see also iGopal-Krishna. Wiita fc Saripallil 
ll989l ; lBlundell fc Rawlingslll999l ). 

At higher redshifts the high flux limit of the 3CRR sam- 
ple leads to an age distribution skewed towards younger ages 
for which the sources are more luminous. This effect plays 
no significant role for the 6CE and 7CRS samples because 
of their lower flux limits. The drop in the source numbers 
in the bin containing the youngest objects is caused by the 
limit we impose on the size of sources to be included in the 
artificial samples. 

The uniform distribution of t in the artificial samples is 
consistent with our assumption of a maximum lifetime com- 
mon to all sources at a given redshift. Any deviations from 
a uniform distribution, apart from those discussed above, 
would imply that our assumption is flawed. However, we can- 
not turn this argument around. The fact that the uniform 
distribution of t 'survives' the source selection in our model 
only demonstrates that, similarly to the situation with the 
distribution of A, the current data does not constrain the 
distribution of source lifetimes. We cannot conclude that 
tmax is the same for all sources at a given redshift. 

We do not initially constrain the distribution of the jet 
power, Qo, to take a specific form as we do with A and t. 
For each source we adjust Qo to give the correct radio lumi- 
nosity using the randomly selected values for A and t. The 
only restrictions we apply are that Qo cannot lie outside the 
ranee from 10 37 W to 5 x 10 40 W. In this way the resulting 
final distribution for Qo arises from the constraints of the 
observed samples themselves. Figure fTUl shows the results for 
the artificial samples related to the 3CRR and 7CRS sam- 
ples for three redshift ranges. We do not consider the result 
for the 6CE sample here as this sample also has an upper 
flux limit which affects the distribution of Qo at the high 
power end. 

The shape of the distribution at the low power end is de- 
termined by the flux limit of the observed samples. Weaker 
jets produce less luminous lobes below the flux limit unless 
the sources are located in very dense environments and/or 
they are young. For high jet powers we expect that virtually 
all objects are above the sample flux limit, independent of 
environment or age, and that the Qo distributions shown 
in Figure \W\ are representative of the entire radio source 
population in the universe. This is supported by the agree- 
ment in the slope of the distributions for the 3CRR and the 
7CRS samples in this region. The distribution for the arti- 
ficial sample corresponding to the 3CRR sample turns over 
for higher jet powers because the flux limit for this sample is 
higher. We do not show results for redshifts beyond z = 0.6, 
because the flux limit influences the high power end of the 
distribution, even for 7CRS. 

The high power part of the Qo distribution is well ap- 
proximated by a power-law. We find that dA/dlogQo oc 
Qo 90 for all three redshift ranges shown in Figure [TO] We 
do not find any evidence for a change of this power-law slope 
as a function of redshift. We have also tested whether the 
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Figure 16. The final distribution of Qo i n three redshift ranges, 
^total is the total source number in each redshift bin. The solid 
line is for 3CRR, the dashed line for 7CRS and the dotted line 
shows a power-law with exponent —1. We have normalized the 
number of the sources so that the curves at the high power end of 
the artificial samples equivalent to 3CRR and 7CRS are aligned. 



slope changes when we change some parameters in the sim- 
ulation proccess. If we adopt (5 — 1.5, we find an exponent 
of —1.05, and the artificial samples from the GLF give an 
exponent of —0.89. Again, if we adopt a different initial dis- 
tribution for A as discussed in Section 8.3. For the Gaussian 
distribution the slope of the Qo distribution becomes slightly 
steeper, i.e. dA/dlogQo oc Qo~ ' 95 ■ For the extreme assump- 
tion of a single value of A for all sources we find a change 
of the power-law exponent to —1.3. We did not attempt a 
formal fit of a power-law to the high power end of the Qo 
distribution as it is not clear below which value of Qo we 
need to exclude data points from the fit to avoid the flux 
limit effects. Given the problems in determining the power- 
law slope, we conclude that dA/dlogQo oc Qq ■ This result 
does not depend on the uncertain value of f3 or the plau- 
sible values for the parameters describing the RLF, or the 
assumed distribution of A as long as very extreme assump- 
tions, e.g. a single value of A for all sources, are avoided. 

Various other authors have tried to constrain the distri- 
bution of Qo as well. Most of these studies present values for 
dN/AQo rather than dA/dlog Qo. We can easily convert our 
result by noting that dTV/dQo = Q ~ 1 dA/dlog Qo oc Qq 2 . 
BRW from their work suggest that the power-law slope 
should be —2.6. More recently BW07 argued for a steeper 
exponent of —3 while KB07 find a value of —1.6. Our result 
is somewhat flatter than those of BRW and BW07 while it 
is slightly steeper than that of KB07. The differences may 
be caused by the steeper luminosity evolution of the BRW 
model and the assumption of a single value for A in BW07. 
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Figure 17. Final distribution of lobe pressures, p c , in two redshift 
ranges. The upper plane is for ^ z ^ 0.2 and the lower plane 
is for 1.0 $C z $C 1.2. iVtotal i s the total source number in each 
redshift bin. The solid line is for 3CRR, the dotted line for 6CE 
and the dashed line for 7CRS. 

KB07 ignored the effect of radiative energy losses of the 
synchrotron emitting electrons on the luminosity evolution, 
which may explain the flatter distribution found by them. 

The KDA model assumes that the energy contents of 
the magnetic field and of the relativistic, synchrotron radi- 
ation emitting particles in the lobe initia lly follow the min- 
imum energy relation (e.g. [Longair 1994). Radiative energy 
losses change this situation somewhat for older sources, but 
the deviations of the model from minimum energy conditions 
are small for most sources. Under minimum energy condi- 
tions the strength of the magnetic field and the volume of 
the emission region completely determine the radio lumi- 
nosity. In our model the volume of the lobe is determined 
by its length and the energy density of the magnetic field 
is proportional to the lobe pressure, p c . Hence the measure- 
ments of lobe length and radio luminosity allow only a small 
range of possible lobe pressures. The data from the observed 
samples should therefore tightly constrain the distribution 
of lobe pressures in our artificial samples. 

Figure [T7] shows the distribution of the lobe pressure 
for two different redshift ranges. The agreement between the 
three samples at low redshifts is good, indicating that the 
different flux limits do not influence the shape of the distri- 
bution. We take this as evidence that at low redshift we ob- 
serve the entire FRII population. At high redshift the sample 
with the highest flux limit, 3CRR, shows a shift of the peak 
in the pressure distribution to higher pressures compared to 
the other samples. Sources with lower lobe pressures are not 
luminous enough to be included in the sample and hence are 
missing. The agreement between the artificial samples corre- 
sponding to 6CE and 7CRS may indicate that these samples 
still include the entire source population at this higher red- 
shift. 

The pressure distribution is peaked and the peak shifts 
to higher values at higher redshifts. Our Model C im- 
plies that the maximum lifetime of a source decreases as 
(1 + z)" 2 ' 4 with increasing redshift. This also implies a pro- 
portional decrease of the average age of sources, (i), included 
in the sample. Our implementation of the KA97 model pre- 



dicts that (p c ) oc (t)~ 2 and therefore (p c ) oc (1 + z) 4 ' 8 . We 
expect the average lobe pressure to increase by a factor of 
roughly 22 (1.35 in the logarithmic scale used in Figure ITT)) 
between redshifts z = 0.1 and z — 1.1. This is consistent 
with the shift of the peak in the pressure distributions for 
the samples with the lower flux limit in the Figure 1171 

Given the good agreement of the pressure distributions 
between the artificial samples corresponding to 6CE and 
7CRS in both redshift ranges, we argue that the samples 
cover the entire FRII source population at both redshifts 
without excluding too many sources through the respective 
flux limits. If so, then the shift of the peak in the pressure 
distribution with redshift is evidence for cosmological evo- 
lution of the population. This conclusion stems from our as- 
sumption of conditions close to those described by minimum 
energy assumptions inside the lobes. It does not depend on 
the additional assumption of Model C of a decreasing max- 
imum lifetime of sources with increasing redshift. The ap- 
plication of Model B would result in the same shift in the 
pressure peak. Unless we invoke a systematic change with 
redshift away from minimum energy conditions inside the 
lobes, the data imply that the average lobe pressure is in- 
creasing rapidly with redshift out to about z = 1. Beyond 
this redshift the flux limits of the samples used here exclude 
some of the FRII-type objects and so we cannot determine 
whether this trend continues. Also, the current data do not 
allow us to determine whether this increase in pressure with 
redshift is due to decreasing source lifetime and/or to an 
increase in the density of the surrounding medium. 

BRW and BW06 use a similar approach to constraining 
the radio source population. We extend the work of BW06 
by testing our model predictions in the complete P-D-z data 
cube. We do not consider the spectral index as the evolution- 
ary models of FRII sources themselves restrict the possible 
range of the spectral index. Our result is similar to that of 
BW06 in that the BRW model gives the worst agreement. 
However, between the other two models, BW06 prefer the 
MK model as it produces a better description of the source 
number ratios at different redshifts. In our work, we use the 
RLF to avoid the fitting of number ratios and the KDA 
model gives a slightly better fit than the MK model. The 
KDA model is also simpler than the MK model and so we 
prefer the KDA model in this paper. Using the RLF, we do 
not need to constrain the distribution of jet powers as BRW 
and BW06 were required to do. The final distribution of Qo 
in our artificial samples agrees with a power-law distribu- 
tion with an exponent of approximately —2, which is flatter 
than the assumption of BRW and the best fit values found 
by BW06. 



11 CONCLUSIONS 

We have constructed a method for generating artificial flux- 
limited radio samples. We use these samples to study the 
cosmological evolution of the FRII source population by 
comparing with observed samples. We use three different 
models for the evolution of the linear size and the radio 
emission of individual FRII sources from KDA, BRW and 
MK. Comparing artificial with observed samples, the 3-D 
KS test indicates that the KDA model gives the best pre- 
diction. The MK model can give acceptable predictions, but 
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it is more complex with more model parameters and it does 
not significantly improve the fitting results. The BRW model 
gives the worst prediction because the radio luminosity de- 
creases too fast as the FRII source get older. 

The properties of FRII sources are required to evolve 
with redshift in order for our artificial sample to fit the ob- 
servations. For /3 = 2, we introduce two models to meet the 
requirements: 

- Model B: The maximum value of the environment 
density parameter A evolves with redshift. We describe 
Model B with best fitting parameters by t — r x 5 x 10 7 yr, 
A = r x 3.7 x 10 17 (1 + z) s ' 8 kgm _1 , where r is a random 
number with uniform distribution between and 1. 

- Model C: The maximum jet age evolves with redshift. 
Similar to Model B, Model C can be expressed by t = r x 
2.7 x 10 7 (1 + z) 2 ' 4 yr, A = r x 6.4 x lQ 17 kgm -1 . 

Both Model B and Model C produce artificial samples in 
agreement with the observed samples according to the 3-D 
KS test. We cannot decide whether a reduced maximum age 
or a denser environment is present at high redshift. Of course 
we cannot rule out a combination of these two effects. In fact, 
the very strong cosmological evolution of the density in the 
source environment required in Model B appears unrealistic. 
It is likely that a more moderate increase of A combines with 
a reduction of the jet lifetimes to produce the observed effect. 

Using our artificial samples, we study the distribution 
of the properties of FRII sources. The distributions of A and 
t are uniform, as we assume initially. The initial setting of 
the distribution 'survives' the source selection in our simula- 
tion and their real distribution cannot be constrained by the 
currently available data. Unlike A and t, we do not initially 
constrain the distribution of Qo, and a power law distribu- 
tion arises naturally from the final artificial sample. We find 
a power law exponent of x ~ —2, and the slope shows no 
significant change at different redshifts up to z = 0.6. We 
also study the distribution of the lobe pressure. The peak 
shifts to a higher value at higher redshift up to z = 1.2. 
This shift arises from our assumption of conditions inside 
the lobes close to those expected from the minimum en- 
ergy requirement. It does not depend on other details of our 
model. 
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